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Abstract 

A  uniform,  dense  film  of  spinel  LiMn204  (0. 1  pm  thick)  has  been  prepared  by  the  electrostatic  spray  deposition  (ESD)  technique.  The 
electroanalytical  behavior  of  this  electrode  is  elucidated  by  application  of  electrochemical  impedance  spectroscopy  (EIS).  The  data  have 
been  modeled  using  an  equivalent  circuit  approach.  An  excellent  fit  was  found  between  measured  data  and  an  equivalent  circuit,  comprising 
Li+  migration  through  surface  film,  potential-dependent  charge  transfer  resistance,  semiinfinite  Warburg-type  element,  reflecting  solid  state 
Li+  ion  diffusion  and  a  finite  space  Warburg-type  element,  describing  both  diffusion  and  accumulation  of  lithium  at  the  very  low  frequency. 
The  apparent  chemical  diffusion  coefficient  of  lithium  in  the  spinel  phase  was  found  within  10  12  <  //  ,  <  10  9  cnr2  s'1  as  a  function  of 
electrode  potential  with  minima  at  the  potentials  corresponding  to  the  voltammetric  peaks.  The  intercalation  capacitance  was  found 
0.7  <  Cl  <  47  mF  cm~2  exhibiting  maxima  at  the  potentials  corresponding  to  the  voltammetric  peaks.  ©  2001  Elsevier  Science  B.V.  All 
rights  reserved. 
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1.  Introduction 

LiMn204  is  an  attractive  cathode  material  for  lithium 
secondary  batteries  [1-6].  Integration  of  such  batteries  into 
electronic  devices  at  the  chip  level  necessitates  fabrication 
of  thin  film  components  of  the  battery.  In  order  to  comply 
with  the  requirements  of  the  microelectronics  fabrication 
technologies,  the  technique  of  preparing  these  films  plays  a 
major  role.  Various  techniques,  such  as  chemical  vapor 
deposition  and  sputtering  have  been  attempted  to  prepare 
thin  films  of  diverse  lithium  battery  components.  Recently, 
Schoonman  et  al.  [7-9]  have  utilized  the  electrostatic  spray 
deposition  (ESD)  technique  to  fabricate  thin  films  of  lithium 
battery  materials  like  LiCoCE,  LiMn204  and  LiNi02.  This 
technique  offers  several  advantages,  such  as  a  simple  and 
low-cost  setup,  low  temperature  synthesis  and  an  easy 
control  of  the  film  composition  and  morphology. 

Recently,  we  have  successfully  prepared  dense  films  of 
spinel  LiMn204  by  ESD  [10].  The  electrochemical  behavior 
of  thus  prepared  films  was  investigated  by  means  of  cyclic 
voltammetry  and  potential  step  chronoamperometry  in  a 
solution  of  1  M  LiC104  in  propylene  carbonate  (PC).  It  has 
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been  found  that  the  oxidation  and  reduction  of  the  film, 
which  are  accompanied  by  the  extraction/insertion  of 
lithium  ions,  occurred  almost  quantitatively  around  4  V 
versus  Li/Li+  with  excellent  reversibility.  The  present  paper 
aims  at  characterizing  the  performance  of  these  electrodes 
by  application  of  electrochemical  impedance  spectroscopy 
(EIS). 

A  major  restriction  of  EIS  characterization  technique  is 
that  multiple  equivalent  circuits  can  be  fitted  to  an  impe¬ 
dance  spectrum  yielding  somewhat  equivocal  results.  None¬ 
theless,  when  the  physicochemical  situation  of  the 
electrochemical  systems  are  clear  and  the  choice  of  the 
model  is  based  on  parallel  morphological  spectroscopic,  and 
electroanalytical  studies,  then  it  may  considerably  reduce 
the  ambiguity  and  enable  a  selected  circuit  analog  to  be 
decided  on.  Another  limitation  concerns  with  the  data 
analysis  procedure  itself.  Complex  linear  least  squares 
(CNLS)  fitting  of  the  data  to  a  theoretical  model  and/or 
equivalent  electrical  circuit  is  the  best  method  of  quantita¬ 
tive  analysis.  Such  fitting  provides  estimates  of  the  para¬ 
meters  and  their  standard  deviations  [11-13].  Unfortunately, 
as  noticed  by  Inzelt  and  Lang  [14],  in  the  majority  of  papers 
no  standard  deviations  or  relative  standard  deviation  of 
the  parameters  are  given.  The  goodness  of  fit  is  illustrated 
in  the  figure  only.  Control  of  the  goodness  of  the  fit  is  also 
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necessary,  i.e.  to  check  whether  the  derived  parameters 
depend  on  the  number  of  elements  or  not,  or  on  the  method 
of  weighting. 


2.  Experimental  details 

LiMn204  cathodes  were  prepared  by  the  electrostatic 
spray  deposition  (ESD).  The  ESD  set-up  used  in  this  work 
was  similar  to  that  described  in  the  literature  [7-9].  Sum¬ 
marily,  a  high  DC  voltage  (in  a  range  of  0-20  kV)  was 
applied  between  an  electrically  conductive  substrate  and  a 
metal  capillary  nozzle,  which  is  connected  to  a  precursor 
solution.  At  an  appropriate  flow  rate,  the  precursor  solution 
was  atomized  at  the  orifice  of  the  nozzle,  generating  the 
spray.  The  spray  moves  toward  the  heated  substrate  under 
electrostatic  force,  and  owing  to  the  pyrolysis  of  precursors, 
a  thin  layer  is  formed  on  the  substrate  surface.  An  ethanol 
precursor  solution  of  25  mM  LiNCE  +  50  mM  Mn(NC>3)2 
was  pumped  at  0.1  ml/h  to  a  stainless  steel  nozzle  (0.8  mm 

i.d.),  which  was  placed  2.5  cm  above  the  gold  flag  substrate, 
heated  typically  at  400°C.  The  spraying  nozzle  was  main¬ 
tained  at  12  kV  against  the  electrode  substrate.  During  the 
film  deposition,  the  surface  of  the  gold  electrode  was  masked 
by  means  of  an  aluminum  foil  in  order  to  limit  the  exposed 
area  to  a  circle  of  6  mm  in  diameter.  The  film  thickness  and 
morphology  were  studied  by  means  of  scanning  electron 
microscope  (JSM-5310LV),  while  the  crystallographic 
structure  was  studied  using  a  Cu  Ka  radiation  with  an 
incident  angle  of  2°  (Shimadzu  XD-D1). 

The  electrolyte  was  1  M  solution  of  LiCICE  in  a  propylene 
carbonate  (PC)  (1  M  =  1  mol/dm3).  The  water  content  of  the 
PC  solution  was  less  than  20  ppm.  Electrochemical  char¬ 
acterization  of  the  cathode  materiel  was  carried  out  using  a 
three-electrode  cell  (T-cells).  Lithium  foils  were  used  as  the 
reference  and  counter  electrodes. 

EIS  experiments  were  performed  using  Solartron’s  1286 
electrochemical  interface  and  125 5 -frequency  response  ana¬ 
lyzer  driven  by  the  Corrware  for  Windows  software  (Scrib¬ 
ner  Associates).  The  impedance  data  generally  covered  the 
frequency  range  of  6  MHz  to  10  kHz  with  an  ac  voltage 
signal  of  ±10  mV.  The  data  were  analyzed  with  a  complex- 
nonlinear-least  square  program  developed  by  Boukamp 
[15].  By  means  of  this  program,  the  optimized  component 
values  for  a  variety  of  different  equivalent-circuit  models 
could  be  rapidly  evaluated.  This  procedure  permits  distin¬ 
guishing  between  alternative  models  and  an  assessment  of 
the  physical  interpretation  of  the  components  used  in  any 
given  model. 

In  the  following,  Zw,  ZFLW  and  ZFSW  are  the  semiinfinite 
Warburg,  finite-length  Warburg  and  finite-space  Warburg 
diffusion  impedance,  respectively.  The  circuit  description 
codes  (CDC)  of  these  elements  are,  W,  O  and  T,  respectively. 

W  =  — (1) 
Yo(Vjcq) 


O  = 


tanh  (Bs/jco) 

Yq(VJm) 


T  =  coth(By//m) 

Yo(V]co) 


(2) 

(3) 


Yq  and  B  are  the  fitted  parameters.  B  depends  on  the 
diffusion  coefficient,  Dl?  and  the  layer  thickness,  /,  it  follows 
that: 


(4) 


The  term  F0  contains  both  the  diffusion  constant,  the  con¬ 
centration,  as  well  as  the  physical  dimension  of  the  system 
under  study. 


3.  Results  and  discussion 

Theoretical  analysis  of  the  complex  impedance  for  charge 
injection  involving  layers  of  finite  thickness  was  developed 
by  Ho  et  al.  [16]  during  the  study  of  lithium  injection  and 
diffusion  in  thin  solid-state  WO3  films.  According  to  the 
authors,  it  is  expected  that  the  impedance  of  the  electroactive 
films  will  exhibit  at  high  frequency  a  charge  transfer  domi¬ 
nated  regime.  At  lower  frequencies,  diffusion  of  charge  in 
the  film  dominates  the  impedance  results.  The  diffusion 
coefficient  and  the  film  thickness  control  the  frequency 
range  of  the  diffusion  behavior.  Finally,  the  finite  thickness 
of  the  film  will  limit  the  extent  of  the  diffusion  behavior  at 
low  frequencies,  and  a  “redox  capacitance”  is  observed. 

Fig.  1  presents  a  family  of  Nyquist  diagrams  and  their 
corresponding  Bode  plots  measured  in  the  range  of  poten¬ 
tials  from  3.80  to  4.40  V  (at  open  circuit  potential  condi¬ 
tions)  along  with  the  anodic  branch  of  the  CV  curve,  i.e. 
during  the  course  of  deintercalation  (due  to  scale  limitation, 
the  low  frequency  at  the  beginning  (3.80  V)  and  the  end 
(4.40  V)  of  the  process  is  not  displayed  in  the  Nyquist  plots). 
If  the  time  constants  for  the  different  processes  are  suffi¬ 
ciently  different,  the  diagram  consists  of: 

1.  A  pronounced  semicircle  at  the  high-to-medium  fre¬ 
quencies. 

2.  In  the  low  frequency  range  (0.79-2  Hz),  a  straight  line 
with  a  slope  of  45°  from  the  real  axis  is  observed  which 
corresponds  to  semiinfinite  Warburg  impedance. 

3.  At  the  very  low  frequencies,  the  45°  begins  to  give  way 
to  a  more  or  less  vertical  line  corresponding  to  finite 
space  diffusion  process  or  the  differential  intercalation 
capacity  of  the  electrode  at  the  measurement’s  potential 
[17,18], 

3.1.  Impedance  analysis:  methodology 

Several  equivalent  circuits  (EqCs)  were  sought  which 
could  describe  the  impedance  response  satisfactorily  under 
all  conditions  examined.  For  the  sake  of  compar  ison,  we  will 
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Fig.  1 .  A  family  of  impedance  plots  of  the  LiMn204  (0. 1  pm  thick  and  0.2  cm2  area)  electrode  for  various  potentials  during  the  course  of  deintercalation  (the 
values  of  the  potentials  are  indicated).  Points  represent  experimental  spectra  and  solid  lines  represent  simulated  curves  according  to  Model  IV  of  Fig.  2.  (a) 
Nyquist  plots;  (b)  admittance  plots  and  (c)  Bode  representation. 


restrict  ourselves  hereafter  in  presenting  only  circuits  that 
are  more  consistent  with  the  data  and  the  physical-chemical 
properties  of  the  system  under  study. 

To  study  the  differences  between  several  possible  models, 
there  are  rigorous  ways  to  get  an  indication  on  how  well  the 
modeling  function  reproduces  the  actual  data  set.  First,  by 
observing  the  parameter  values  and  their  relative  error 
estimate  (%).  The  Chi-squared  function  (y2)  also  provides 
a  good  indication  of  the  quality  of  the  fit.  Generally,  the  y1 
statistic  should  be  less  than  10  6  if  the  data  are  said  to  fit 
within  the  noise  level  of  the  measurement.  Therefore,  a 
value  of  10~5  ~  10~4  indicates  a  reasonable  good  fit.  y1 
based  on  modulus  weight  factor  [19]  is  defined  here  as 


follows: 

N  (  yf  yf  .  /  yii  7//  \2 

2  V-"1  vG  G,sim/  '  i  ^fsim / 

X  “  4^  Z[ 2  +  Z"2  (  } 

where  Z[  and  Z"  are  real  and  imaginary  parts  of  experimental 
data,  Z[sim  and  Z"sim  are  real  and  imaginary  parts  of  simu¬ 
lated  data. 

Then,  the  actual  data  set  can  be  observed  in  a  graph  of  the 
relative  residuals,  A Z[  and  A Z"  versus  log  (o.  Following 
Boukamp  [19]  the  residuals  are  defined  by: 
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where  co  —  27i/is  the  angular  frequency  of  the  ac -perturba¬ 
tion,  IZsiml  is  the  vector  length  of  the  modeling  function. 

3.2.  Expected  models 

In  the  absence  of  any  surface  layers,  the  impedance  of  an 
insertion  electrode  can  be  interpreted  in  terms  of  a  simple 
network  consisting  of  the  bulk  electrolyte  resistance,  Re,  in 
series  with  a  RWC  network  made  of  two  parallel  current 
branches:  (i)  a  faradaic-current  branch  represented  by  a 
charge  transfer  resistance,  Rct,  serially  connected  to  a  War¬ 
burg  impedance  Zw,  and  (ii)  a  nonfaradaic-branch  consisting 
of  the  double  layer  capacitance,  Cd|.  This  is  the  well-known 
Randles-circuit.  An  interface  with  a  surface  layer  can  be 
modeled  by  adding  another  RslCsl  component  (where  “si” 
stands  for  surface  layer)  to  the  original  RctCdi  network  [11] 
(Fig.  2  Model  I).  Previously,  Model  I  in  Fig.  2  has  been 
applied  for  example  to  lithium  insertion  in  Sn-oxide  film,  in 
Ti-oxide  film,  in  transparent  Ti-Ce  oxide  films  [20-22]  and 
into  LiMn204  and  cobalt-doped  LiMn2Od  [23].  In  our  work, 
the  use  of  constant  phase  element  (CPE)  has  been  circum¬ 
vented  deliberately  because  our  spectra  were  well  modeled 
by  ideal  capacitors.  Therefore,  correlation  to  practical  sys¬ 
tems  was  preserved  with  practical  physical  significance. 
Model  III  in  Fig.  2  is  basically  similar  to  the  one  applied 
by  Aurbach  et  al.  [24,25]  to  Li  intercalation  process  into 
graphite  and  into  transition  metal  oxides  as  LitCo02, 
LivNi02  and  LivMn204  spinel  electrodes.  ZFLW  is  a  finite  - 
length  Warburg  type  element  reflecting  the  Li  solid-state 
diffusion.  The  Cint  is  a  capacitor  describing  the  accumula¬ 
tion  of  the  intercalant  (lithium)  into  the  host  materials.  We 
have  to  emphasize  that  in  Aurbach ’s  model  there  are  four 
serially  connected  RIIC  networks  describing  a  multisurface 


layer.  Their  circuit  is  known  as  a  combination  of  the  Voigt- 
type  analog  and  the  generalized  Lrumkin  and  Melik  Gay- 
kazan  (LMG)  impedance.  On  the  other  hand.  Model  III  in 
Lig.  2  applies  for  a  one  surface  layer.  We  have  modified 
Model  I  by  incorporating  a  finite  space  Warburg-type  ele¬ 
ment  ZFSW  (Model  IV  in  Lig.  2).  The  diffusion  element, 
ZFSW,  describes  the  diffusion  in  a  medium  where  the  inter¬ 
face  precludes  the  flow  of  the  species  and  is  represented  by 
hyperbolic  cotangent  function  [15,26,27], 

3.3.  Modeling  results  and  interpretation 

Table  1  provides  a  comparison  of  modeling  quality  of  the 
different  models  proposed  in  Lig.  2  when  applied  to  impe¬ 
dance  spectra  obtained  at  4.02  V.  It  is  plainly  seen  from  the 
y2-values  that  Model  I  do  not  describe  adequately  the 
system.  Models  II-IV  seem  to  refine  to  low  /2-values  with 
reasonable  values  for  the  physical  processes  represented. 
However,  a  visual  inspection  of  the  graphical  output  pro¬ 
duced  by  the  fitting  program  (Lig.  3)  shows  that  only  Models 
II  and  IV  were  capable  of  mimicking  the  spectra.  It  is  also 
seen  from  Lig.  4  that  error  distributions  versus  frequency  are 
the  lowest  (<2%)  for  Models  II  and  IV. 

We  have  also  fitted  our  data  to  the  complete  Voigt-type 
analog  and  the  generalized  Lrumkin  and  Melik  Gaykazan 
(LMG)  impedance  used  by  Aurbach  et  al.  [24,25] .  The  fitting 
results  gave  a  good  y2  of  2.10  5.  However,  the  error 
estimates  for  the  RIIC  components  were  considerably  high 
(10-20%).  Moreover,  we  found  no  congruous  variation  of 
the  circuits  components  with  the  applied  potential.  This 
might  be  due  to  the  closeness  of  the  time  constants  of  the 
circuits,  and/or  some  of  the  components  have  only  a  small 
contribution  to  the  overall  frequency  dispersion. 
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Fig.  2.  Equivalent  circuit  representations  of  LiMmCC  film  electrode  considered  in  this  paper  (see  details  in  the  text). 


Table  1 


Fitting  parameters  for  models  represented  in  Fig.  2  applied  to  impedance  diagram  recorded  at  4.02  Va 


RJQ. 

Cdl/mF 

rjq 

Csl/mF 

RJQ 

ZWY0  x  10~3/ 
Q"1  s0-5 

Zflw  To  x  10  3/ 

Q  1  s0-5 

^FLw/s0'5 

-Zfsw  Yq  x  10  3/ 
Q-1  s0'5 

^FSw/s<)'5 

Qnt/mF 

x2 

Model  I 

106  (10.8) 

4.5  (18.5) 

361  (8.09) 

3.13  (136) 

21.2  (85) 

2.39  (8.1) 

- 

- 

- 

- 

- 

4.29  x  10~2 

Model  II 

108  (0.81) 

4.9  (0.86) 

426  (0.86) 

4.44  (7.5) 

35.1  (7.8) 

5.63  (2.6) 

- 

- 

- 

- 

7.7  (1.34) 

4.32  x  10~4 

Model  III 

108  (1.18) 

4.9  (3.31) 

433  (1.29) 

5.1  (10.9) 

36.4  (12) 

- 

6.98  (4.26) 

2.6  (6.28) 

- 

- 

6.3  (1.49) 

9.75  x  10~4 

Model  IV 

107  (0.61) 

4.8  (1.45) 

412  (0.66) 

4.5  (6.35) 

29.8  (5.8) 

6.5  (3.05) 

- 

- 

7.04  (6.3) 

2.12  (5.8) 

- 

2.04  x  10^4 

a  Numbers  between  brackets  represent  error  estimates  (%)  for  the  individual  components  of  the  equivalent  circuit.  /iFI  W  and  5fsw  are  defined  by  Eq.  (4)  in  the  text. 
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Fig.  3.  Measured  impedance  spectra  for  the  LiMmC^  film  electrode  obtained  at  4.02  V  compared  to  spectra  calculated  following  the  models  shown  in  Fig.  2. 
(a)  Nyquist  plots;  (b)  admittance  plots  and  (c)  Bode  representation. 
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Fig.  3.  ( Continued ). 


Phase  angle/Degree 


100 


M.  Mohamedi  et  cil.  /  Journal  of  Power  Sources  93  (2001)  93-103 


Fig.  4.  Relative  residual  plots  for  the  fits  presented  in  Fig.  3.  Real  error 
and  imaginary  error  were  calculated  according  to  Eq.  (6). 

A  tentative  interpretation  of  Model  IV  can  be  made  by 
comparison  with  Goodenough  et  al.  [11]  analysis.  Accord¬ 
ing  to  the  authors,  two  different  physical  processes  may  be 


involved  to  account  for  the  addition  of  an  id  1C  to  the  Randles 
equivalent  circuit:  (i)  adsorption  of  Li+  or  propylene  car¬ 
bonate  onto  the  surface  of  the  electrode  without  charge 
transfer  and  (ii)  formation  of  an  ionically  conducting  but 
electronically  insulating  surface  layer  at  the  electrode  sur¬ 
face.  To  distinguish  between  these  two  models  Goodenough 
et  al.  [11]  have  studied  the  time  dependence  of  the  ac- 
impedance  response  and  the  variation  of  the  circuit  para¬ 
meters  with  applied  voltage.  Their  results  have  shown  that 
only  the  surface  layer  model  contains  equivalent  circuit 
parameters  that  vary  in  a  self-consistent  manner  with  the 
electrochemical  processes  they  represents  as  the  cell  con¬ 
ditions  are  varied. 

Conclusively,  in  agreement  with  Goodenough  et  al.  [11], 
we  ascribe  the  RS\CS\  to  the  resistance  for  Li+  ion  migration 
through  the  surface  film  and  film  capacitance,  as  probably 
formed  by  electrolyte  decomposition  product.  It  is  not  trivial 
then  to  assign  the  second  RWC  to  the  electron  transfer 
resistance  of  the  intercalation  reaction  and  capacitance  of 
a  double  layer.  Although  Model  II  describes  adequately  the 
observed  response,  the  fitting  results  were  further  evaluated 
for  all  applied  potentials  according  to  Model  IV  using  ZFSW 
instead  of  Cmi  because  of  the  straightforward  physical 
significance  of  the  former  parameter.  Results  of  fitting  are 
exhibited  in  Fig.  1. 

3.4.  Model  parameters 

We  observe  that  Rct  is  high  for  high  potentials  (3.8  and 
4.4  V)  (Fig.  5).  The  charge  transfer  resistance  being  high  for 
high  potentials  is  consistent  with  the  observations  on  the 
voltammograms  where  the  current  is  zero  (see  Fig.  6,  here- 


E/V  vs.  Li/Li+ 

Fig.  5.  Plots  of  the  charge  transfer  resistance,  Rct,  and  the  surface  layer 
resistance,  Rs\  vs.  potential.  Rcl  and  Rs\  were  obtained  by  fitting  the 
experimental  impedance  spectra  with  Model  IV  of  Fig.  2. 
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Voltammogram  with  0.5  mV  s1 
•  I=C  *v 

L 


Fig.  6.  Comparison  of  the  voltammetric  response  of  a  LiMn204  electrode 
at  0.5  nrV/s  (solid  line)  with  ClXv  (points).  Cl  is  the  film  redox 
capacitance  obtained  from  Eq.  (14)  and  v  is  the  scan  rate  used  in  obtaining 
the  actual  voltammogram. 


after).  The  surface  layer  resistance,  Rs\  was  found  more  or 
less  constant  in  the  4.02-4. 15  V  potential  range  (Fig.  5).  The 
values  of  Rsi  at  4.40  and  3.80  V  have  been  left  since  the  error 
estimates  obtained  in  the  fitting  on  the  Rs\  and  Csi  parameters 
were  extremely  large. 

The  values  of  the  surface-layer  capacitance  were  in  the 
range  of  23-32  mF  cm-2  (Fig.  7),  in  agreement  with  those 
reported  by  for  the  pure  spinel  [23] .  The  C&\  values  ranging 
from  23-34  mF  cm”2  (Fig.  7)  are  characteristics  of  double 
layer  capacitance  and  support  the  equivalent  circuit  inter¬ 
pretation. 

3.5.  Determination  of  the  exchange  current  density 


The  charge  transfer  resistance  is  related  to  the  exchange 
current  density,  i0,  by  the  equation: 


Ret  — 


RT 

nFio 


(7) 


The  exchange  current  densities  were  calculated  using  Eq.  (7) 
for  various  applied  potentials.  The  results  are  shown  in 
Fig.  8.  The  exchanged  current  density  lies  between  0.01 
and  0.4  mA/cm2.  We  found  that  i0  exhibits  a  maximum  at  the 
potential  of  the  fust  CV  peak  and  then  monotonously 
decreases  for  higher  potentials. 


3.6.  Determination  of  Li-ion  chemical  diffusion  coefficient 


E/V  vs.  Li/Li+ 

Fig.  7.  Plots  of  the  double  layer  capacitance,  C<ji,  and  the  surface  layer 
capacitance,  Csi  vs.  potential.  Cji  and  Csi  were  obtained  by  fitting  the 
experimental  impedance  spectra  with  Model  IV  of  Fig.  2. 


•  Under  semiinfinite  conditions  (eo  >->-  2D\A/12)  in  the  case 
when  the  time  constants  of  the  different  physical  pro¬ 
cesses  occurring  at  the  electrode  are  separated,  Du  can  be 
obtained  from  an  analysis  of  the  Warburg  impedance,  Zw 
[16],  which  is  expressed  in  the  complex  plane  as: 


Zw  = 


V„ 


d  E 


sflnFAD d* 


(1  -j)co  1/2  =  <7(1 


-  Oco”1/2 


(8) 


where  <x  is  the  Warburg  coefficient,  Vm  the  molar  volume  of 
the  host  structure,  d E/Ax  the  slope  of  the  coulometric  titra- 


E/V  vs.  Li/Li+ 


The  chemical  diffusion  coefficient  Du  may  be  obtained 
by  either  several  ways: 


Fig.  8.  Exchange  current  density,  iQ  vs.  potential;  iQ  was  determined  by 
means  of  Eq.  (7). 
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tion  curve  versus  mobile  ion  concentration  x  at  each  x  value, 
/  the  maximum  length  of  the  diffusion  pathway,  and  n  the 
number  of  equivalent  (n  —  1).  The  Warburg  coefficient  can 
be  obtained  from  the  slope  of  the  Z"  versus  co-1/2  or  Z! 
versus  eo_1/2. 

In  the  general  case,  the  Warburg  impedance  in  the  faradaic 
branch  must  be  related  to  the  diffusion  coefficients  of  the 
insertion  ion  in  both  the  electrolyte  and  the  electrode.  In  our 
experimental  spectra,  the  observed  semiinfinite  Warburg 
impedance  might  be  confounded  with  diffusion-controlled 
process  in  solution  and  in  the  oxide.  In  order  to  remove  all 
doubts  about  the  meaning  of  the  Warburg  behavior,  we  have 
studied  the  influence  of  electrolyte  concentration,  especially 
in  the  low  frequency  range  related  to  lithium  diffusion  in  the 
oxide.  No  significant  variation  in  the  Warburg  slope  with  the 
electrolyte  concentration  was  observed  in  the  range  of  1  M 
to  10  mM  LiC104.  This  indicates  that,  for  our  experimental 
conditions  (1  M  LiCICU),  lithium  diffusion  in  the  bathing 
electrolyte  is  not  a  limiting  factor,  and  proves  that  the 
semiinfinite  Warburg  behavior  observed  in  our  experimental 
spectra  is  relevant  to  lithium  diffusion  in  the  solid  phase. 

•  If  the  homogeneous  diffusion  (co  yy  2Du/l 2)  is 
observed  (last  part  of  the  diagram),  the  current  is  approxi¬ 
mately  90°  out  of  phase  with  voltage.  Under  this  condi¬ 
tion,  the  real  part  of  the  impedance  is  independent  of  the 
frequency.  The  extrapolated  real  intercept  of  the  straight- 
line  region  gives  rise  to  the  limiting  low  frequency 
resistance,  RL,  defined  by  Ho  et  al.  [16]: 


Z' 


ym  d  E  l 
nFA  ck  3 Du 


(9) 


The  authors  also  express  the  imaginary  part  of  the 
impedance  as: 


1  _  Tm  d  E  1 

coCl  nFA  dv  col 


(10) 


1  _  d Z" 
Cl  dco-1 


(ID 


E/V  vs.  Li/Li+ 

Fig.  9.  Semi-logarithmic  representation  of  the  variation  of  the  chemical 
diffusion  coefficient  of  Li-ions  vs.  potential  for  the  LiMrbCh  electrode. 


The  determination  of  RL  is  not  precise  and  an  inaccuracy 
is  introduced  in  the  calculation  of  the  diffusion  coefficient 
using  Eq.  (12).  More  accurate  values  of  the  Du  can  be 
deduced  from  either  the  semiinfinite  diffusion  (Eq.  (8))  or 
the  finite  space  diffusion  domain  using  Eq.  (13). 

Since  the  Warburg-type  impedance  was  assessed  as  being 
film  diffusion  process  where  the  diffusion  layer  thickness,  /, 
can  replace  the  thickness  of  film  (0. 1  mm)  through  which  the 
diffusion  takes  place.  Therefore,  the  ratio  If  \f  Z\,  in  Eq.  (13) 
was  determined  from  the  fitted  parameter  B  (Eq.  (4)).  Fig.  9 
shows  that  10  “12  <  Du  <  10  9  during  deintercalation 
process.  As  also  demonstrated  in  this  figure,  log  Du  versus 
E  exhibits  two  distinct  minima  at  the  anodic  CV  peak 
potentials. 

3.7.  Determination  of  the  intercalation  capacitance,  CL 


Thus,  in  the  range  of  the  impedance  phase  angle 
approaches  nil,  the  value  of  the  redox  capacitance,  CL 
can  be  determined  graphically  from  the  slope  of  Eq.  (11). 
Therefore,  Du  may  be  deduced  from  Eqs.  (9)  and  (11) 


Du  = 


l 2 

3  RlCl 


(12) 


•  From  the  so-called  finite  space  diffusion  Warburg  impe¬ 
dance  [15,16,26-28]  as  this  element  reflects  both  diffu¬ 
sion  and  accumulation  (capacitive  behavior  at  the  very 
low  frequencies): 


Zfsw  — 


-1/2 

coth 


-1/2 


(13) 


The  low  frequency  capacitance,  was  estimated  using: 


(14) 


Direct  comparison  of  the  intercalation  capacitance,  CL 
determined  from  EIS  data  with  the  voltammetric  response 
can  be  made  by  multiplying  CL  by  the  scan  rate,  v  employed 
in  the  voltammetric  experiments.  In  theory,  these  data  are 
exactly  equivalent  to  that  obtained  from  CV  response  of  an 
ideal  film  electrode  at  sufficiently  slow  scan  rates  (condi¬ 
tions  where  thin  layer  electrochemical  behavior  is  observed) 
[29],  or  that  the  film  be  sufficiently  thin  to  insure  that  the 
effects  of  diffusional  charge  transport  in  the  film  phase  is 
absent.  In  our  previous  study  [10],  we  found  that  that  current 
peak  was  strictly  proportional  to  scan  rate  only  at 
v  <  0.5  mV/s.  Fig.  6  demonstrates  good  agreement  of  1-E 


where  t  =  l2 /Du  is  the  diffusion  constant. 
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curve  reconstructed  from  CL  data  with  the  corresponding 
cyclic  voltammogram  obtained  at  a  0.5  mV/s  scan  rate. 

4.  Conclusions 

Characterization  of  LiMn204  thin  film  made  by  ESD 
technique  has  been  achieved  by  means  of  impedance  spec¬ 
troscopy.  Experimental  spectra  have  been  modeled  follow¬ 
ing  rigorous  criteria.  The  derived  model  contains  a  RIIC 
network  related  to  surface  layer  via  which  Li+  migrates,  a 
Randles-circuit  describing  the  electron  transfer  and  the 
lithium  ions  solid  state  diffusion  and  a  finite  space  diffusion 
element  describing  both  lithium  diffusion  and  accumulation. 
Specific  parameters  of  interest  determined  in  this  investiga¬ 
tion  of  LiMn204  thin  films  were:  (i)  diffusion  coefficient,  (ii) 
redox  capacitance,  (iii)  surface  film  capacitance,  resistance, 
and  (iv)  exchange-current  density.  An  effective  method  of 
testing  robustness  of  the  model  will  be  to  change  the 
experimental  conditions  (electrolyte,  temperature,  film 
thickness,  etc.)  to  observe  how  well  the  model  predicts 
the  changes,  i.e.  in  order  to  find  physically  reasonable 
dependencies  that  help  to  discriminate  adequately  between 
the  derived  parameters.  This  will  be  the  topic  of  our  next 
paper. 
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